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!.  (U)  Introduction 

The  purpose  of  this  work  was  to  assess  experimentally  the  feasibility  of  locating  power  generators  by  time- 
correlating  their  acoustic  energy  relayed  by  three  or  more  spatially-separated  acoustic  sensors.  Correlation  tech¬ 
niques  arc  very  commonly  used  in  the  process  of  localizing  noise  sources^1'11*  and  many  of  the  factors  affecting 
sound-ranging  arc  well  understood112*. 

The  specific  problem  addressed  in  this  work  was  eliminating  the  spurious  maxima  from  the  correlation 
function  of  the  generator  signal.  These  spurious  maxima  result  from  the  dominant  periodic  components  existing  in 
the  generator’s  acoustic  signal.  Depending  on  the  nature  of  the  sensor  configuration  and  the  particular  generator 
tested  ten  or  more  maxima  of  equal  amplitude  may  appear  in  the  correlation  function.  This  makes  it  difficult  to 
correctly  identify  the  true  correlation  maximum. 

This  report  discusses  a  successful  technique  to  eliminate  the  effect  of  periodic  components  from  the  correlation 
function  and  defines  relatively  simple  methods  of  correctly  identifying  these  types  of  targets.  Chapters  II,  III,  IV, 
and  V,  respectively,  di'cuss  the  theory  behind  this  experiment,  the  experimental  apparatus  and  microphone  con¬ 
figuration,  the  experimental  analysis  and  results,  and  summarize  the  findings/ 


II.  (U)  Theory 


The  location  of  sound  sources  is  done  by  using  time-of-arrival  techniques.  To  establish  the  coordinates 
of  the  acoustic  source  by  this  method  the  following  must  be  known:  (a)  a  minimum  of  two  independent  delr,  times 
(time  differences  between  the  same  event  received  at  several  microphone  locations);  and  (b)  the  corresponding 
microphone  separating  distances  (base  lines).  In  practice,  the  microphone  array  configuration  is  known  and  the  only 
parameters  to  be  measured  are  the  delay  times  between  the  elements  of  at  least  a  threc-nicrophonc  array.  If  the 
acoustic  source  is  of  impulsive  nature  the  time  difference  between  signals  received  at  several  microphone  locations 
is  easily  defined.  This,  however,  is  not  so  in  the  case  of  continuous  signals.  In  the  latter  case  cross-correlation  tech¬ 
niques  can  achieve  results  analogous  to  those  obtained  for  the  pulse  case. 

A.  Correlation  Analysis 

The  cross-correlation  of  two  waveforms  is  a  measure  of  their  similarity  as  a  function  of  the  time  shift  between 
them.  Mathematically  this  can  be  achieved  by  computing  the  cross-correlation  function.  That  is,  the  two  signals  are 
multiplied  and  summed  over  a  specified  time  window  T.  This  is  done  for  different  values  of  time  r.  For  the  value  of 
r  at  which  the  two  signals  match,  a  maximum  value  of  the  correlation  function  is  obtained.  This  value  of  r  is  equal 
to  the  delay  time  of  one  signal  relative  to  the  other.  The  mathematical  expression  for  the  discrete  cross-correlation 
function  is: 


Rij  (m,k)  =  dj 


rsrm 


Xj(n)  X:  (n-k) 


Here,  Xj(n)  is  the  signal  on  channel  i  at  the  n^  sample  instant,  Xj(n-k)  is  the  signal  on  channel  j  at  the  (n-k)th 
sample  instant,  k  is  the  delay  between  channels  in  samples,  N  is  the  correlation  window  in  samples,  and  m  is  the 
sample  instant  for  start  of  the  correlation  window. 


In  the  case  of  discrete  plotting  of  a  series  of  measurements  in  the  form  of  pairs  of  observation  points  Xj(n'» 
and  Xj(n),  a  normalized  measure  of  correlation  lying  between  -1  and  +1  is  often  desirable.  The  digital  form  of  ,'~r 
normalized  correlation  is:  _ r  _ _  — - - J  „ 

r  N+m  1  r  N+m — 1N2 - 


Pij(mW=Rij(m,k).  UX)X2(n)  -  I  £ 


The  normalized  correlation  function  was  suggested  as  a  method  to  reduce  spurious  correlation  maxima. 
rationale  behind  it  is  as  follows:  Write  Eq  (1)  as, 

Rjy  (m;k)  =  Oj(m)  Oj(m-k)  pjj(m,k) 
f  N+m  i  1/2 


N+m 

«!<">>-  i;Sxi 

n=m 


If  there  is  a  large  change  in  energy  in  the  signal  from  one  correlation  record  to  the  next,  then  it  is  possible 


Rij  (mj  Jcj)  >  Rjj  (m2.  k2) 


(5) 


while, 


kl)  <  Pij(m2>  k2> 


(6) 


simply  because 

0j(mJ)0j(,T'i  -kj)  »£7i(m2)aj(m2-k2)  (7) 

If  the  signal  is  quasi-stationary  in  the  sense  that  its  spectral  shape  does  not  change  with  time,  then  py  will  be 
independent  of  the  record  time  m.  Hence  it  will  be  a  more  consistent  estimator  of  the  true  correlation  maxima  than 
fly  if  the  signal  energy,  above,  is  non-stationary. 

For  an  example,  consider  the  application  of  correlation  techniques  to  an  acoustic  receiving  system  consisting 
of  a  three-element  airay.  The  coordinates  of  the  array  relative  to  a  known  origin,  in  this  case  a  control  processing 
center,  will  be  known. 

B ,  Time  of  Arrival  Techniques 

In  Figure  1  the  field  deployment  of  such  an  array  is  illustrated.  It  has  been  mentioned  above  that  the  signal 
from  the  source,  whose  position  relative  to  the  array  is  S  (X,Y),  arrives  at  different  times  at  Ml,  M2,  and  M3.  The 
difference  in  time-of-arrival  (delay  time  r)  can  be  determined  using  cross-correlation  techniques.  Using  the  delay 
time  to  locate  the  source  of  disturbance  one  must  find  a  relationship  between  r  and  the  geometry  of  the  system. 

To  simplify  the  calculations  we  have  assumed  a  planar  geometry  (i.e.,Z=Zj  =  Z2  =  Z3).  If  this  is  not  the  case  an 
array  of  at  least  four  sensors  is  required  to  locate  the  sound  source. 


Figure  1.  Location  of  a  Sound  Source  by  Triangulation 
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Using  an  orthogonal  coordinate  system  the  microphone  configuration  can  be  designated  as  M 1  (X j ,  Yj), 
M2(X2,  Y2),  and  M3(X3,  Y3).  Using  simple  geometrical  configurations  the  following  relationships  can  be  obtained 
for  the  time  delays  of  the  signal  produced  by  an  acoustic  source  S  having  coordinates  (X,  Y)  and  recorded  at  M 1 , 
M2,  and  M3: 


atl*»2*  £  ^(YrY)2  +  (X-X1)2J  1/2 - 

rCY2-Y)2  +  (X2-X)2l  1/2j, 

[[(vWx-xJ2]  >'2- 

£(Y- Y3)2  +  (X3-X)2  11 
*3  ~  *2  “  y  [[(Y-Y3)2  +  (X3-X)2]  1/2- 
J(Y2-Y)2  +  (X2-X)2j  1/2J  . 


In  the  above  equation  V  is  the  velocity  of  sound  in  the  atmosphere  and  is  given  by: 


■  Vt  * v" 


Here  V0  =  332  meters  per  second  at  0  °C,  T  is  the  absolute  temperature  i.i  degrees  kelvin,  T0  is  the  ambient 
temperature  in  degrees  kelvin  and  Vw  is  the  velocity  of  the  wind. 

- Equations  8  are  the  equations  of  hyperbolas  whose  foci  are  the  fixed  points  MI7  M2,  and  M 3.  Therefore, 

for  a  given  t  we  can  determine  a  hyperbola,  that  is  a  set  of  points  whose  difference  in  distance  from  the  fixed 
points,  for  example,  Ml  and  M2,  is  equal  to  rj2  (a  constant).  Two  such  hyperbolas  should  intersect  at  least  once. 
The  point  of  intersection  determines  the  location  of  the  source  of  disturbance.  As  the  source  changes  position 
the  delay  times  will  also  change  and  hence  a  new  intersection  will  be  defined.  Therefore,  by  computing  the  delay' 
times  between  different  microphone  pairs  in  real  time  the  source  of  disturbance  can  be  tracked. 

C.  Signal  Prewhitening 

Signal  prewhitening  is  used  primarily  to  reduce  the  periodic  structure  of  the  correlation  function.*10*  In  the 
present  analysis  the  prewhitening  and  cross-ccrrelation  scheme  consists  of  taking  two  time  functions  Xj(t)  and 
X2(t),  taking  the  discrete  Fourier  transforms  Fj(w)  and  F2(w),  multiplying  the  first  with  the  conjugate  of  the 
second  in  order  to  find  the  cross-spectral  power  density,  prewhitening  the  result  in  order  to  suppress  periodicity, 
and  then  taking  the  inverse  transform  to  obtain  the  cross-correlation.  Note  that  if  the  prewhitening  step  is  omitted 
the  result  is  a  normal  cross-correlation  function.*13* 

For  example,  if  Xj(t)  and  X2(t)  are  two  time  functions  their  Fourier  transforms  Fj(w)  and  F2(«)  can  be 
written: 

FiM^a^wJ  +  jbjM 

and 

F2(w)=a2(w)  +  jb2(w)  (10) 
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Then  the  cross-spectral  power  density  is: 

S  [ ,  2(w)  =  F  «(co)  F  *2(w)  *  [a  i  (w)  ai  (w)  +  b  |  (w)b2(oj)] 

+  j  (a2(w)  b  1  (co)  -  b2(w)a  j  (w)] . 

Any  prewhitening  scheme  in  the  frequency  domain  can  be  considered  to  be  a  multi'  lication  of  S  j  2(u>)  by  a  function 
K-(cj)/10)  In  general,  there  are  no  restrictions  placed  on  k2(w).  However,  in  order  to  preserve  phase  information, 
K-(co)  is  required  to  be  positive  real.  After  performing  the  multiplication  we  have, 

(co)  . 

Ti>2M  =  S1i2k2(co)  (11) 


K2(W)  = 


2[a2+02j  1/2 


a  *  1  j  (to)  32  (w)  +  b  j  (co)  b2  (co)  and 

0  =  a2  (to)  b  j  (co)  -  b2  (co)  a  j  (co)  (13) 

At  this  point  take  the  inverse  transform  of  Tj  2(00).  Using  the  convolution  theorem  the  inverse  transform  of  Tj  2(«) 
equals  the  inverse  transform  of  Sj  2(00)  convolved  with  the  inverse  transform  of  K2(co).  Then  since  the  inverse 
transform  of  S]  2(co)  is  the  cross-correlation  Rj  2(f)  of  X](t)  and  X2O),  the  result  is  the  corss-correlation  of 
Xj(t)  and  X2(t)’ convolved  with  a  function  k(r),’where  k(r)  is  the  inverse  transform  of  K2(co). 

D.  Spectral  Analysis 


Spectral  analysis  is  based  upon  computation  of  the  Fourier  integral: 


/"> 


X(t)  exp  (-2  n  ft)  dt 


In  order  to  understand  how  computation  of  the  Fourier  integral  helps  to  make  meaningful  measurements,  consider 
the  Fourier  transform  written  in  its  sine-cosine  form: 


sX(0  «/° 


X(t)  (cos  2ffft  -  j  sin  2rrft)  dt 


This  equation  states  that  the  transform  averages  a  time  function  input  X(t)  with  a  set  of  sines-cosines  to  determine 
the  content  of  X(t)  at  some  frequency  f.  Thus  the  transform  resolves  the  time  function  into  a  set  of  components  at 
various  frequencies  much  as  a  set  of  analog  filters  would.  However,  it  not  only  yields  the  amplitude  of  each  frequency, 
but  also  resolves  the  in-phase  (real,  cosine)  component  and  the  quadrature  (imaginary,  sine)  component,  thereby 
giving  magnitude  and  phase  information  which  is  difficult  to  obtain  in  any  other  way. 
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III.  (U)  Experimental  Apparatus  and  Experimental  Arrangement  (U) 

A .  (C)Target  -  Sensor  Configuration  (U) 

(C)  The  target,  50  Hz,  75  KVA  generator  motor  set,  was  tested  in  a  semi-jungle  environment  area  on  Auxiliary 
Field  No.  6  at  Eglin  AFB,  Florida.  The  test  was  conducted  in  order  to  establish  criteria  for  locating  and  identifying 
generator  type  targets.  Figure  2  depicts  the  relationship  between  the  sound  source  and  the  sensors.  The  detection 
system  consisted  of  modified  Phase  III  commikes  (MCOM),  hand-emplaced  commikes  (HEC),  and  an  RADC  acoustic 
system.  The  Phase  III  commikes  were  connected  to  a  central  recording  station  via  an  RF  link,  while  the  RADC 
microphones  were  hard-wired  into  a  recording  system. 

(C)  The  MCOM  and  HEC  were  placed  at  points  designated  by  stake  markers  on  the  four  legs  (north,  south,  east, 
and  west)  from  the  source.  The  stakes  were  in  100-meter  increments  up  to  1  kilometer  and  from  there  in  500-meter 
increments  up  to  2  kilometers.  The  source  as  indicated  in  Figure  2  was  positioned  off-center.  The  location  was 
approximately  20  meters  from  center  on  the  north  leg  and  approximately  8  meters  from  center  on  the  west  leg.  It 
should  be  noted  that  the  exhaust  was  directed  parallel  to  the  east  leg. 

(C)  The  two  MCOMS  and  one  HEC  were  collocated  on  each  of  the  four  legs.  Two  MCOMS  allowed  for  backup 
data  in  case  a  MCOM  failed  during  a  test.  The  HEC  gave  continuous  audio  data  for  evaluation  of  the  audio  signal 
at  any  time  during  the  operation  of  the  MCOMS. 

(C)  The  RADC  acoustic  microphones  were  positioned  at  various  points  along  the  four  legs  of  the  test  array. 

Figure  3  depicts  the  four  different  arrays  used  by  RADC  and  the  positions  of  the  microphones  relative  to  the  source. 

B. (C)  Sensor  System  Description  (Phase  III  Commikes)  (U) 

(C)  The  HEC  transmitted  continuous  audio  while  the  MCOM  transmitted  type  I  messages  upon  recognition  of  the 
source.  The  MCOM  v/as  commanded  to  the  audio  mode  at  various  times  throughout  the  test.  The  MCOM  consists  of 
a  Phase  Ul  commike  modified  by  NADC  (Naval  Air  Development  Center).  The  commike  was  modified  to  consist 
of  detec  tion  circuitry  which  was  included  in  the  audio  module,*14)  as  shown  in  Figure  4. 

(C)  The  detection  circuitry  consisted  of  a  phase  lock  loop  which  sweeps  from  72  to  78  Hz  nominal.  A  signal 
within  this  frequency  range  must  be  present  at  least  10  seconds  (enable  timer)  before  the  unit  will  begin  to  process 
the  signal.  After  the  10-second  period,  the  signal  must  be  present  for  7  minutes  more  (verification  timer),  with  no 
dropouts  exceeding  10  seconds.  At  the  end  of  7  minutes  a  type  I  message  is  transmitted  signifying  the  presences  of 
the  “source.”  If  the  signal  should  dropout  for  more  than  10  seconds  the  entire  process  starts  over  with  the  enable 
timer.  At  any  time  during  the  processing  of  the  signal  the  MCOM  can  be  commanded  into  the  audio  mode. 

(C)  There  is  a  certain  finite  amount  of  time  during  which  any  MCOM  can  transmit  audio  after  being  commanded 
into  the  audio  mode.  The  amount  of  continuous  audio  time  after  a  command  to  the  audio  mode  is  20  seconds 
nominal.  Since  the  unit’s  battery  life  is  limited  mainly  by  the  number  of  times  it  transmits  audio,  the  audio  command 
is  used  infrequently. 

(C)  The  operating  characteristics  of  the  MCOM  are  the  same  as  the  Phase  III  commike  but  will  be  restated  here 
for  clarity.  The  frequency  response  of  the  MCOM  is  from  50  to  2  KHz  with  output  distortion  less  than  10%  for  the 
audio  mode.  The  output  noise  level  is  24  millivolts  rms  maximum  while  the  instantaneous  dynamic  range  is  32  db. 
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BOTE:  Exhaust  from  source 
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parallel  to  east  leg 
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Figure  2.  Source  and  Sensor  Locations  at  Auxiliary  Field  No.  6  Eglin  AFB  FL 
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Figure  4.  Phase  III  Commike  Modification  for  PAVE  ONYX  Mission 
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(C)  The  information  transmitted  by  the  MCOM  is  received,  digitized,  and  transmitted  (S-band)  by  an  airborne 
segment  data  link  to  a  receiving  and  analyris  center.  The  digitizing  rate  is  approximately  4  KHz.  The  digital 
information  is  then  transmitted  using  an  S-band  transmission.  The  link  is  depicted  in  Figure  S.  The  audio  digital 
data  is  then  transformed  back  to  analog  data  (D/A)  at  the  receiving  and  analysis  center. 


Figure  5.  S-Band  Data  Link  for  Phase  III  Commikes 

C.  (U)  RAPC  Acoustic  System 

I  Tne  RADC  experimental  apparatus  consisted  of  the  sound  system,  microphones,  the  micrometeorological 

instrumentation,  ar.d  the  recording  equipment. 

f  1.  Sound  System  —  The  sound  system  was  a  100-watt  system  consisting  of  two  loudspeakers  and  an  audio 

tape  recorder.  The  speakers  were  Lansing  Models  LE  1 54  which  have  a  relatively  flat  response  in  the  frequency 
range  40-6000  Hz.  The  audio  portable  recorder  was  a  Tandberg  Model  1 1  with  frequency  response  40-20,000  Hz. 
In  addition  to  the  generator  target  the  following  sound  sources  were  used  in  this  experiment: 

a.  Broad  Band  Noise  -  Filtered  noise  in  octave  bands  were  used  as  a  noise  source.  The  bands  are  40-80 
Hz,  80-160, 160-320, 320-640, 640-1280  Hz.  The  duration  of  these  signals  is  15  seconds  each.  These  signals  were 
used  to  determine  the  ability  to  cross-correlate  on  pure  noise  and  to  compare  these  results  with  cross-correlation 

|  functions  obtained  from  the  “source  truck’’  signals.  In  addition  to  noise-like  signals,  truck  signatures  contain  strong 

spectral  lines.  v  .  ; 

b.  Swept  Tones  —  These  are  swept  tones  and  pure  tones  covering  the  frequency  range  between  50-150  riz. 

i  The  structure  of  these  signals  is  as  follows:  (1)  50  Hz  for  15  seconds;  (2)  50-150  Hz  range  swept  in  15  seconds;  (3) 

1 50-50  Hz  swept  in  1 5  seconds;  (4)  50  Hz  pure  tones  for  1 5  seconds.  These  signals  were  included  in  order  to  make 
direct  comparison  with' the  source  signals  which  appear  to  sweep  the  same  frequency  range.  Also,  it  was  of  interest 
I  In  itself  to  test  cross-correlation  on  these  types  of  signals. 
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c.  Moving  Target  (Truck)  -  The  signatures  from  a  moving  truck  along  a  prescribed  path  have  been 
recorded.  These  signatures  were  used  as  an  interfering  target. 

d.  Interfering  Targets 

(1)  In  this  case  three  trucks  were  used  at  the  same  time.  The  trucks  were  separated  by  100m. 

(2)  Prop  and  jet  aircraft  were  used  as  interfering  targets  in  various  configurations.  The  masking  experi¬ 
ments  are  designed  to  determine  the  effectiveness  of  cross-correlation  techniques  to  locate  and  follow  targets  of 
interest  as  a  function  of  base  line  in  the  presence  of  interference  phenomena. 

2.  Microphones  -  In  answering  the  question  concerning  the  applicability  of  cross-correlation  techniques  for 
a  stationary  periodic  target,  four  arrays  of  microphones  were  used.  These  are  depicted  in  Figure  3.  An  array  con¬ 
sisted  of  four  omnidirectional  microphones.  The  microphones  are  GR-P1 560-953 l’s  and  have  a  frequency  response 
which  is  flat  in  the  range  between  2  and  6000  Hz. 

3.  Micrometeorological  Instrumentation  -  A  meteorological  system  manufactured  by  Meteorology  Research, 
Inc.  (MRl)  was  used  in  this  experiment.  The  instrumentation  was  mounted  on  a  6-foot  tower.  The  system  is  capable 
of  giving  the  temperature,  wind  speed,  and  wind  direction.  Using  these  data,  information  about  temperature  stability 
and  the  turbulence  in  the  area  can  be  obtained. 

4.  Recording  Equipment  -  The  signals  from  the  microphones  are  hardwired  to  a  central  recording  system. 

The  lines  are  driven  by  a  set  of  amplifiers  which  are  gain  switchable.  The  total  dynamic  range  of  the  system  is  90  db 
and  the  gain  can  be  switched  between  27  db  and  97  db.  The  steps  are  10  db  each.  The  bandwidth  of  the  system  is 

2.  to  2000  Hz.  The  output  of  the  amplifiers  was  fed  into  a  14-ciiannel  analog  tape  recorder.  The  recorder  was  a 
Honeywell  5600  magnetic  tape  recorder.  Out  of  the  14  channel"  in  the  FM  record  mode  four  were  used  for  the 
microphones  and  three  were  used  for  the  meteorological  channels.  The  edge  track  A  was  used  for  voice  commentary. 
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IV.  (U)  Analysis  and  Results 


The  experimental  data  for  this  effort  were  collected  at  Eglio  AFB  Florida.  The  experimental  details  and  the 
type  of  the  sound  source  {referred  to  as  “source”  in  the  following  sections)  against  which  the  data  were  collected 
are  summarized  in  Chapter  III  of  this  report.  The  analysis  of  the  data  was  performed  by  using  the  following  equip¬ 
ment:  (a)  Federal  Scientific  Correlation  and  Probability  Analyzer; (b)  Federal  Scientific  UA-6B  Spectrum  Analyzer; 
and  (c)  PDP-9  General  Purpose  Digital  Computer. 

A.  Spectrum  Analysis 

To  determine  the  properties  of  the  “source,"  its  sound  intensity  and  the  frequency  spectrum  of  the  emitted 
signal  were  examined. 

1 .  Sound  Intensity  -  To  establish  a  safe  distance  over  which  the  “source”  could  be  detected  the  sound- 
pressure  level  (SPL)  one  meter  away  from  the  target  was  examined.  The  SPL  at  any  distance  from  the  sound  source 
can  be  determined  using  the  following  relationship: (16),  (I7) 

1(R)  =  (Ws  +  D)  -  20  log  R  -  MR  -  T  (dB)  (16) 

Here,  Ws  is  the  source  acoustic  power  in  dB  relative  to  10'^  watts;  D  is  the  directivity  index;  M  is  the  attenuation 
coefficient  in  dB  per  unit  distance;  and  T  is  the  transmission  loss  in  dB.  The  term  20  log  Represents  the  geometrical 
spreading  loss. 

The  “source”  exhibited  some  degree  of  directionality.  The  SPL  one  meter  away  and  in  the  direction  of  the 
exhaust  system,  or  the  two  sides  adjacent  to  the  exhaust,  was  found  to  be  108  dB,  while  the  SPL  on  the  side  opposite 
to  the  exhaust  was  about  100  dB.  The  background  noise  at  the  test  site  was  of  the  order  of  56  dB  during  the  day  and 
about  40  dB  during  the  night.  The  Eglin  terrain  was  heavily  wooded  and  a  loss  coefficient  O7)  M  =  .0264  dB/in  is 
estimated  for  this  terrain  and  for  this  source  (i.e.,  most  of  the  energy  of  the  source  is  found  to  bo  below  150  Hz). 

- -  — -  =  Using  Equation  (16),  the  SPL  of  the  “source,”  and  the  environmental  noise,  an  effective  range  for  a  signals _ 

to-noise  ratio  equal  to  unity  is:  (a)  in  the  direction  of  the  exhaust:  (1)  daytime,  R=340m  and  (2)  nighttime,  R=460m; 
(b)  in  a  direction  opposite  to  the  exhaust:  (1)  daytime,  R=230m  and  (2)  nighttime,  R=3S0m.  These  results  are  in 
excellent  agreement  with  the  experimentally  observed  values.  However,  data  processing  improved  the  above  detection 
ranges  considerably  (i.e.,  a  factor  of  two  ur  more  is  a  realistic  range  with  minimum  processing). 

t  The  above  analysis  shows  that  a  safe  range  for  detecting  the  source  of  interest,  with  minimum  data 

!  processing,  is  about  600m.  Larger  ranges  could  be  achieved  if  the  background  noise  is  reduced  and/or  by  additional 

'  processing  of  the  data.  The  periodic  nature  of  the  source  of  interest  (see  next  section)  suggests  that  ranges  of  several 

J  kilometers  could  be  achieved  by  averaging  the  data  over  longer  periods  of  time  (i.e.,  several  minutes). 

2.  Source  Spectrum  -  Using  Equation  (14)  the  spectrum  of  the  source  signal  has  been  calculated.  In  Figure 
'  6.  a.  a  typical  power  spectrum  over  an  analysis  range  0-2000  Hz  is  shown.  In  Figure  6.b.  the  source  spectrum  is 

shown  over  the  analysis  range  of  100-2000  Hz  after  the  signal  was  amplified  by  20  dB.  This  chart  was  produced  by 
passing  the  signal  through  a  high  pass  filter  and  then  amplifying  it  by  20  dB.  Figure  6.c.  shows  the  power  spectral 
’■  density  of  the  background  noise  (generator  off),  produced  under  the  same  conditions  as  Figure  6.b.  In  addition  to 

:  the  fact  that  the  source  is  highly  periodic,  Figuie  6  indicates  that  most  of  the  acoustic  energy  is  found  in  the  0-100 

Hz  frequency  band.  Figure  7  shows  details  of  the  generator  spectrum  in  this  frequency  range.  It  is  observed  that  the 
|  signal  has  a  fundamental  of  12.5  Hz  and  that  the  strongest  frequency  component  is  at  25  or  37.5  Hz  (i.e.,  the 

second  or  the  third  harmonic).  Whether  the  25  Hz  or  the  37.5  Hz  line  is  the  strongest  line  in  the  spectrum  depends 
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Figure  6.  The  Power  Spectrum  of:  (a)  the  generator  in  the  analysis  range  0-2000  Hz,  0  dB  gain;  (b) 
the  generator  in  the  analysis  range  100-2000  Hz,  20  dB  gain;  and  (c)  the  background  noise  (generator 
off)  in  the  analysis  range  0-2000  Hz,  20  dB  gain.  (Array  No.  2) 


Figure  7.  The  Power  Spectrum  of  the  Generator  for  Different  Microphones  in  the  Analysis  Range  of  0-100  Hz; 
(a)  in  the  direction  of  the  exhaust,  exciter-on;  (b)  in  the  direction  of  the  exhaust,  exciter-off;  (c)  behind  the 
exhaust,  exciter-on;  (d)  behind  the  exhaust,  exciter-off. 


on  the  direction  of  the  exhaust  pipe  relative  to  the  detector.  That  is,  if  the  micropnone  is  in  the  same  direction  as 
the  exhaust  pipe  outlet,  the  37.5  Hz  line  is  the  strongest  line  in  the  received  spectrum,  while  if  the  microphone  is 
in  the  opposite  direction  to  the  exhaust  pipe  outlet,  the  25  Hz  line  is  the  strongest. 

From  the  mechanics  of  the  source,  the  12.5  Hz  line  appears  to  be  associated  with  the  cylinder  firing  rate. 
The  cylinder  firing  rate  (CFR)  for  a  six-cylinder,  four-cycle  engine  is  given  by:(18) 


CFR=  ^^1  _  [Revolutions Per  Minute]  ^zj 

2x60  [60  -^1x2 

L  minj 

The  engine  tested  had  an  RPM=1500.  From  Equation  (17)  one  sees  that  the  CFR=12.5  Hz.  The  12.5  Hz  line  appears 
to  be  modulating  all  its  harmonics.  This  modulation  is  clearly  seen  in  the  time  function  Figure  8,  where  the  two 
adjacent  peaks  are  80  milliseconds  apart.  This  is  a  distinct  characteristic  of  this  type  of  sound  source  and  can  be  used 
to  correctly  identify  it.  Less  than  one  second  of  data  is  required  to  correctly  classify  the  target.  Notice,  however, 
that  the  12.5  Hz  modulation  is  reduced  when  the  source  is  running  with  no  load  on  (i.e.,  when  the  exciter  is  turned 
off).  A  direct  comparison  of  the  time  function  of  the  source  with  the  exciter  on  and  the  exciter  off  is  shown  in 
Figure  9.  Figure  10  shows  a  typical  time  series  spectrum  of  the  source  with  the  exciter  on  and  the  exciter  off.  When 
the  exciter  is  off  the  lines  are  shifted  slightly  into  higher  frequency  (see  also  Figure  7).  The  behavior  of  the  source 
with  the  exciter  on/off  is  not  well  understood.  Puzzling  also  is  the  fact  that  the  25  or  37.5  Hz  line  is  more  intense 
than  the  engine  firing  rate  (EFR),  which  in  this  case  is  given  by:(18) 


EFR  =  (RPM  x  (6  cyQ)  a  1 S00  x  6  _  75  Hz 

[6oi2£.lx  2  60x  2  (sec) 

I  mini 

This  may  be  due  to  the  mechanics  of  the  muffler.  No  data  have  been  collected  with  the  muffler  off  using  this 
experimental  apparatus  to  verify  this  assumption. 


Another  important  feature  of  the  source  is  shown  in  Figure  6.  In  this  figure,  although  there  is  significant 
energy  in  the  frequency  domain  between  100-2000  Hz  this  energy  is  about  20  dB  less  than  the  energy  found 
between  0-100  Hz.  However,  randomness  in  the  source  signal  is  more  ptonounced  in  the  analysis  range  100-2000  Hz. 


From  the  above  it  is  concluded  that  the  class  of  targets  can  be  easily  identified  on  the  basis  of  the  modula¬ 
tion  associated  with  the  cylinder  firing  rate  and/or  using  the  rich  harmonic  content  and  periodic  nature  of  the 
source  spectrum.  Also  the  fact  that  there  is  enough  randomness  in  the  source  signal  allows  it  to  be  used  as  a  basis 
for  correlation.  In  light  of  this  information  an  attempt  was  made  to  locate  the  source  using  time-of-arrival 
information. 

B.  Correlation  Analysis 

Using  time-of-arrival  techniques  to  locate  a  source  of  disturbance  is  briefly  discussed  in  the  theory  section  of 
this  report.  It  can  be  shown  that  if  an  array  of  three  or  more  sensors  is  receiving  signals  from  a  sound  source 
simultaneously,  the  essential  parameter  to  measure  is  the  difference  between  times-of-arrival  of  the  signal  at  the 
various  receiver  locations  (delay  times).  Knowing  the  difference  between  the  times-of-arrival  of  a  signal  at  several 
locations  and  the  distance  between  the  receivers  (base  line),  the  source  of  disturbance  can  be  located  from  simple 
geometrical  considerations  (see  Section  II. B.).  The  most  common  method  for  determining  delay  times  is  a  technique 
known  as  cross-correlation.  The  necessary  prerequisite  for  this  method  to  work  is  that  the  signal  is  not  periodic.  If 
the  signal  is  periodic  (see  theory)  the  solution  is  not  unique  and  the  technique  completely  fails  unless  time  differences 
less  than  the  period  of  the  signal  are  expected.  This  implies  sensor  separation  less  than  one  wavelength  long.  However, 
this  restriction  is  operationally  difficult  to  realize  with  the  present  IGLOO  WHITE  Phase  III  sensors.  Operationally, 
for  the  present  application,  one  works  with  sensor  separation  which  is  much  larger  than  the  wavelength  of  the 
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AMPLITUDE 


Figure  8.  Time  Waveform  Generated  by  the  Source  as  Recorded  at  the  Four  Microphones  of  Array  Number  2 
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Figure  9.  Time  Waveform  Generated  by  the  Source:  (a)  exciter-on;  (b)  exciter-off 

minimum  frequency  component  in  the  “source”  signal,  it  is  therefore  required  that  the  “source”signals  must 
contain  some  degree  of  randomness  for  this  technique  to  work  at  all. 

f  -Result's  and  Interpretation . .  '  . — 

The  correlation  analysis  summarized  here  is  performed  using  array  number  4.  Calculations  from  the  other 
arrays  are  omitted  since  they  give  similar  results.  The  different  array  configurations  are  discussed  in  Chapter  III. 
Figures  11,  12,  13,  14,  15  and  16  show  the  correlation  function  that  has  been  computed  using  the  output  from 
different  microphones  of  array  configuration  number  4.  In  these  figures  a  display  of  the  cross-correlation  function 
of  the  source  (generator),  the  band-limited  noise,  and  the  swept  tones  are  shown  for  comparison  purposes. 

Table  I  summarizes  the  results  of  this  correlation  analysis.  The  discrepancies  observed  between  theoretical 
and  experimental  values  of  the  delay  times  are  believed  to  be  due  to  the  fact  that  the  theoretical  values  were 
computed  on  the  assumption  that  the  terrain  was  flat.  In  actuality  this  was  not  the  case.  On  the  other  hand,  the 
differences  observed  between  the  calculated  values  of  the  delay  times  of  the  source  and  the  test  signals  using  the 
same  microphone  pair  are  attributed  to  the  fact  that  the  test  signals  when  played  back  into  the  system  were  physically 
separated  from  the  source  by  about  three  meters.  That  is,  a  maximum  delay  difference  in  the  order  of  8  milliseconds 
is  possible.  This  is  well  within  the  differences  observed.  This  point  clearly  shows  the  resolution  of  the  technique. 

The  good  agreement  between  theory  and  experiment  demonstrates  that  a  source  of  disturbance  can  be  localized 
with  a  high  degree  of  accuracy  using  correlation  techniques  applied  over  long  base  lines/11)  The  computations 
summarized  here  were  performed  using  a  Federal  Scientific  Correlation  and  Probability  Analyzer. 

From  the  above  it  should  be  noted  that  although  in  the  analysis  range  between  zero  and  100  milliseconds 
no  real  difficulties  existed  in  reducing  excessive  periodic  components  from  the  correlation  function,  in  all  cases  it 
was  possible  to  produce  the  true  correlation  maximum.  The  real  problem  for  larger  delays  is  by  far  more  complex. 
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Figure  10.  Time  Series  Spectrum  of  the  Source  Signal:  (a)  exciter-on;  (b)  exciter-off 


Figure  12.  Cross-correlation  Function  of  Ml  vs  M3:  (a)  source  signal;  (b)  40-80  Hz  band-limited  white  noise; 

and  (c)  40-150  Hz  swept  tones.  ti 
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1 5.  Cross-correlation  Function  of  M2  vs  M4:  (a)  source  signals;  (b)  40-80  Hz  band-limited  noise;  and 

(c)  40-1  SO  Hz  swept  tones. 


TABLE  I.  Correlation  Analysis  (Array  No.  4) 


EXPERIMENTAL  DELAY 

SOURCE  NOISE 
r  (msec)  r  (msec) 

500  500 

528  535 

90  96 

28  36 

412  408 

438  442 

When  the  analysis  range  is  extended  beyond  the  1 00  milliseconds  range,  additional  spurious  peaks  made  it  very 
difficult  to  locate  the  true  maximum.  Figures  1 1,  12, 13, 15,  and  16  clearly  demonstrate  this  point.  In  these 
figures  the  correlation  functions  of  the  source  and  different  test  signals  using  different  microphone  pairs  are 
summarised.  Notice  that  using  the  signal  as  it  is,  the  additional  maxima  occurring  every  80  milliseconds  arc  very 
difficult  to  eliminate.  These  periodic  maxima  are  due  to  the  12.5  Hz  signal  modulation  whose  nature  is  discussed 
in  an  earlier  section  of  this  report. 

Faced  with  this  dilemma  two  possible  solutions  were  discussed:  (a)  deploy  a  sensor  system  whose  base  lines 
would  be  less  than  30  meters;  or  (b)  find  the  origin  of  the  periodic  components  and  determine  means  to  reduce 
their  effect  on  the  correlation  function.  Since  the  first  solution  was  very  difficult  to  realize  from  the  operational 
point,  of  view,  emphasis  was  placed  in  the  area  of  signal  processing  with  the  possibility  of  extracting  features 
(parameters)  that  were  characteristic  of  the  random  components  in  the  signal.  This  was  hopefully  based  on  the 
assumption  that  by  its  nature  any  mechanical  system  produces  a  signal  which  in  addition  to  the  periodic  components 
contains  also  some  degree  of  randomness.  That  is, 

f(t)=*  ap(t)+  bn(t)  (19) 


MICROPHONE  PAIR 

Ml  vs  M2 
Ml  vs  M3 
M 1  vs  M4 
M2  vs  M3 
M2  vs  M4 
M3  vs  M4 


THEORETICAL  DELAY 
r  (msec) 

510 

536 

78 

26 

432 

458 


where  the  ap(t)  and  bn(t)  are  the  periodic  and  random  parts  of  the  emitted  signal  fft).  Based  on  this  assumption 
the  correlation  function  becomes, 


ll ,2 <T>  "  ^  jf1  f  1  (0  f2(l  - t)  dt  *  yj  T a,  a2  p !  (t). 
T 

p2(t-r)dt+—  f  bjb;n(t)n(t-r)dt 
•'-T 


The  question  that  arises  now  is.  what  is  the  contribution  of  the  second  term  to  the  correlation  function?  In 
answering  this  question  the  following  analysis  was  performed.  Over  a  fixed  window  the  auto-correlation  function 
of  the  source  signal  was  computed  after  the  signal  was  passed  through  a  1 00-1000  Hz  high  pass  filter.  The  back¬ 
ground  signal  was  also  auto-correlatcd  under  identical  conditions.  The  results  of  this  analysis  are  summarized  in 
Figure  17.  By  comparing  the  two  auto-correlation  functions  (source  on/source  off)  it  can  be  shown  that  the 
additional  energy  at  r°0  In  the  case  cf  the  source  suto-corrclation  function  is  due  to  the  contribution  of  the 
random  components  In  the  source  signal.  Notice  that  the  100-1000  Hz  filter  was  used  to  eliminate  the  contribution 


Figure  17.  Auto-correlation  of  100-1000  Hz  band-limited:  (a)  source  signal;  (b)  background  noise. 

of  the  strong  periodic  components  of  the  signal  that  are  found  below  100  Hz  (Figures  1 1  and  12).  The  effect  of 
the  filtering  on  the  cross-correlation  function  is  shown  in  Figure  18.  For  this  analysis  the  signals  of  microphones 
M3  and  M4  have  been  used.  Notice  that  for  the  case  where  the  1 50-1000  Hz  filter  was  used  the  spurious  maxima 

have  been  greatly  reduced  and  a  discernible  true  correlation  maximum  is  observed  at  r=438  milliseconds.  This - 

figure  clearly  demonstrates  that  signal  preconditioning  is  required  before  cross-correlating  signals  produced  by 
generators  and  other  generator-like  targets.*10) 

In  addition  to  the  above,  the  following  analysis  was  performed  with  the  help  of  a  general  purpose  digital 
computer:  (a)  the  computation  of  the  normalized  cross-correlation  function,  and  (b)  the  computation  of  the  cross¬ 
correlation  function  after  prewhitening  the  signals.  Prewhitening  is  used  here  with  the  intent  to  reduce  the  effect 
of  the  periodic  components  of  the  signal  in  ihe  correlation  function.*101  The  net  effect  would  be  similar  to  tint 
demonstrated  in  Figure  18  where  special  filtering  techniques  were  applied. 

Using  Equation  (2)  and  a  digital  computer  the  normalized  cross-correlation  function  has  been  computed. 
The  results  of  this  analysis  are  shown  in  Figure  19.  In  the  majority  of  the  runs,  although  the  normalized  correlation 
function  was  highly  periodic,  the  true  correlation  maximum  was  the  strongest  maximum  in  the  function,  as  was 
expected.  For  reasons,  however,  that  we  cannot  explain,  In  some  correlation  rum  this  was  not  the  case.  In  the  case 
of  Figure  19  the  signals  of  microphones  M7,  and  M3  have  been  cross-correlated  producing  a  maximum  at  r*28 
milliseconds.  From  array  number  4  and  Table  I  we  see  that  this  maximum  corresponds  to  the  true  time  difference 
between  times-of-arrival  of  the  source  signal  at  microphones  M2  and  M3.  In  most  of  the  runs,  similar  results  were 
obtained  between  the  other  microphone  pairs.  However,  since  this  analysis  did  not  give  reproducible  results  100% 
of  the  time,  it  was  decided  to  apply  prewhitening  techniques  before  the  signals  were  cross-correlated.  The  theory 
behind  the  prewhitening  scheme  is  summarized  in  Chapter  11.  The  results  of  this  analysis  are  shown  in  Figure  20. 
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Figure  19.  Normalized  Cross-correlation  M2  vs  M2  of  Unwhitened  Source  Signals  (r=28  msec). 


Figure  1 


In  this  figure  the  signals  of  microphones  M2  and  M3  used  to  produce  Figure  14  were  prewhitened  and  then  were 
cross-correlated.  In  this  case,  as  in  the  one  above,  the  maximum  at  r=28  milliseconds  is  the  true  maximum  and 
stands  well  above  the  spurious  maxima.  The  analysis  done  with  other  pairs  of  n.icrophones  using  prewhitening 
techniques  as  a  function  of  integration  window  is  shown  in  Figure  21 .  In  this  figure  the  correlation  maximum  is 
centered  over  the  analysis  window.  Using  this  technique  all  computer  runs  produced  discernible  correlation  maxima. 
On  the  other  hand  notice  that  the  periodic  component  has  been  greatly  reduced  as  compared  to  Figure  10. 

Once  a  workable  correlation  technique  was  established  an  attempt  was  made  to  locate  the  source  using 
data  that  were  obtained  through  an  operational  system  (see  Chapter  111).  Although  the  IGLOO  WHITE  Phase  III 
system  filtered  out  most  of  the  energy  below  50  Hz  because  of  a  notch  filter  used  at  about  60  Hz  (Figure  22),  the 
correlation  maximum  could  be  found  quite  easily  (Figure  23).  Using  the  delay  times  (r)  between  different  micro¬ 
phone  pairs  and  Equation  8  one  can  solve  for  the  coordinates  of  the  source.  These  results  indicate  that  mechanical 
sources  whose  main  signal  is  highly  periodic  could  be  located  using  correlation  techniques. 

Notice  that  in  addition  to  the  method  described  in  this  report  for  identifying  the  generator-type  targets, 
a  technique  based  strictly  on  looking  for  the  engine  firing  rate  has  been  successfully  deployed  in  SEA.U9)  ^  brief 
description  of  this  type  of  detection  is  given  in  Chapter  HI  of  this  report.  For  more  details  see  reference  ( i  9). 
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Figure  22.  Time  Series  Spectrum  of  the  Source  Signal  Using  Data  Obtained  from  an  Operational  System 

(Eglin  AFB). 
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V.  (C)  Summary  and  Conclusions  (U) 


(C)  This  work  was  concerned  with  application  of  acoustic  base  line  techniques  to  detection,  identification,  and 
location  of  enemy  weapon  systems  powered  by  generators.  An  acoustic  base  line  system  locates  by  measuring  time- 
of-arrival  of  acoustic  signals  produced  by  the  target  of  interest,  in  this  case  the  generator,  at  three  or  more  acoustic 
sensors  distributed  in  a  relatively  broad  array.  Such  time-of-arrival  (TOA)  techniques  have  been  used  successfully  in 
other  applications.  In  this  case,  however,  the  long  base  line  (acoustic  sensor  separation)  and  the  highly  periodic  nature 
of  acoustic  signals  generated  by  the  source  lead  to  severe  ambiguity  in  source  location.  The  technical  challenge  in  this 
work  involved  separation  of  the  periodic  components  of  the  source  signals  and  preservation  of  the  more  random 
components  for  use  in  arrival  time  measurement  and  subsequent  target  location. 

(C)  Experimental  work  was  performed  at  Eglin  AFB,  Florida  and  was  designed  to  provide  information  on  all 
variables  of  interest  including  weather  and  interference.  The  sources  used  were: 

1.  Foreign  generators 

2.  Pure  and  swept  tones 

3.  Band-limited  white  noise 

4.  Explosions 

5.  Aircraft  interference 

6.  Vehicles  singly  and  in  groups. 

(C)  Detection  was  performed  using  the  following: 

— - 1.  Hard-wired  microphones  wired  into  an  RADC  recording  facility. . *•—*—**.* . 

2.  Existing  IGLOO  WHITE  system  with  Phase  III  modified  commike  collocated  with  RADC  microphones. 

All  aspects  of  the  experiment  including  design,  installation,  and  operation  of  the  instrumentation  system  were 
performed  by  RADC  personnel  with  support  from  3246TW,  Eglin  AFB,  Florida. 

(C)  Subsequent  to  the  field  phase  of  the  program,  the  data  were  analyzed  at  RADC  using  correlators,  spectrum 
analyzers,  and  digital  data  processing  equipment.  The  results  of  the  analysis  successfully  demonstrated  the 
following: 

1.  Detection  and  location  of  the  source  was  possible  over  base  lines  of  600  meters  by  applying: 

a.  Signal  prewhitening  techniques;  and/or 

b.  Stochastic  correlation  techniques. 

2.  Accuracies  within  a  few  meters  were  achieved  under  calm  weather  conditions. 

3.  Identification  of  the  target  was  performed  based  on  the  modulation  of  the  acoustic  signal  by  the  fundamental 
frequency  of  the  source. 
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(C)  The  technique  described  above  provides  an  approach  to  remote  monitoring  of  enemy  lines  of  communications 
and  fixed  installations  for  tactical  target  interdiction  and  strategic  purposes.  The  most  important  current  example  is 
location  of  SAM  generators  in  a  jungle  environment.  Presently  this  type  of  source  is  very  difficult  to  locate  using 
other  detection  techniques.  This  work  not  only  provides  a  tool  for  accurately  locating  SAM  generators  with 
available  hardware,  it  also  provides  the  means  for  properly  classifying  the  target. 
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APPENDIX 


The  time  domain  correlation  program  was  written  in  FORTRAN  IV  for  use  on  a  PDP-9  general  purpose 
digital  computer  by  RADC/DCT1.  The  correlation  function  is  normalized  (Eq  2)  and  printed  out  in  integer  format. 
A  plot  of  the  normalized  correlation  function  can  be  obtained  by  a  software  plotting  program. 

The  correlation  program  has  the  capability  of  handling  835  points  per  correlation.  The  resolution  depends  on 
the  digitizing  rate  of  the  data  in  conjunction  with  the  835  point  per  correlation. 

The  three-sensor  target  location  program  was  also  written  in  FORTRAN  IV  by  RADC/DCTI.  The  array 
geometry  and  the  time  delays  of  the  acoustic  signatures  of  interest  between  the  three  sensors  are  needed  for  a 
location. 

The  general  equations  used  in  the  digital  program  for  the  three-sensor  target  location  using  TOA  are: 

V0 r  1 2  *  (Xj  +  Y J)  1 11  -  !(X0  -  X2)2  +  (Y0  -  Y2)2]  1  /2 

0  0  (21) 

v0r13  =  (^  +  Y2)  I/2-[(Xo-X3)2  +  (Y0.Y3)2J  1/2 

°  (22) 

where 

Vq  -  speed  of  sound. 

r  12  -  Delay  time  for  the  acoustic  signature  of  interest  between  sensor  1  and  sensor  2. 
r  |3  -  Delay  time  for  the  acoustic  signature  of  interest  between  sensor  1  and  sensor  3. 

Xq,  Yq  -  Cartesian  coordinates  of  the  target  locatic  l. 

X2,  Y2  -  Cartesian  coordinates  for  sensor  2. 

X3,  Y3  -  Cartesian  coordinates  for  sensor  3. 

The  location  of  sensor  1  is  at  the  origin  of  the  Cartesian  coordinate  system,  therefore  Xj  =  0  and  Y  j  *  0. 

Solving  for  the  target  location,  the  following  equations  are  obtained: 

Xq  *  A  +  BR0  (23) 

Y0  a  A,  +  Bj  Rq  (24) 


Ro-(X2  +  Y2)  1/2 

A  *  Y3  (X2  +  Y2)-Y2(X2  +  Y2)-Y3(V0r12)2  +  Y2(V0r13)2 
2(X2Y3-X3  Y2) 

B"Y3(VoT12)-Y2(V0r]3) 

X2Y3-X3Y2 


AI 


(28) 


A,  =  X2(X2  *  Y|)  *  X3(X2  +  Y2)  +  X3(V0  t12)2  -X2(V0r13)2 
2(X2  Y3-X3Y2) 


i 


Bi3X2(V0t13)-X3(V0t12) 

x2y3-x3y2 


(29) 


Squaring  equations  (23)  and  (24),  then  adding  the  two  equations  together,  the  following  quadratic  solution  for  Rq 
is  formed: 

Rq  =  -J±  yJ2-4»K 


2N 


where 


J  =  2(AB  +  Aj  Bj) 
N  =  B2  +  Bj2-1 
K  -  A2  +  A}2 


(30) 

(31) 

(32) 

(33) 


As  long  as  the  target  is  within  the  triangle  formed  by  the  three  sensors,  a  unique  solution  is  obtained  from  equation 
(30).  R0  is  then  substituted  into  equations  (23)  and  (24)  for  the  Cartesian  coordinate  of  the  target  location  relative 
to  sensor  1. 


The  block  diagram  of  Figure  24  depicts  the  s*eps  involved  when  signal  prewhitening  is  used  to  reduce  the 
periodic  structure  of  the  correlation  function.  The  tiieory  and  equations  used  to  develop  the  digital  computer 
programs  for  signal  prewhitening  are  described  and  presented  in  Section  H.C.  Below  is  a  description  of  the  symbols 
in  the  block  diagram  of  Figure  24. 

—  X]  (t),  X2  (t)  —  Time  varying  functions  _  - 

Fj  (to),  F2  (to)  -  Discrete  Fourier  transform  of  Xj  (t)  and  X2  (t)  respectively. 

F2*  (to)  -  Complex  conjugate  of  F2  (to). 

K2  (to)  -  Prewhitening  functions 
FFT  —  Fast  Fourier  Transform 
p  j  2(t)  -  Cross-correlation  function 

WSF  -  Window  Shaping  Function  (i.e.  Hanning,  Hamming  etc.). 


M 


O  U  o  o 


C  THREE  SENSOR  TARGET  LOCATION  PROGRAM. 

C  ROME  AIR  DEVELOPMENT  CENTER. 

C  LT  MARK  D.  DICKINSON  (DCTI). 

C 

C 

C  CALL  INIT  INITIALIZES  THE  TEKTRONIX  6ll  DISPLAY. 

CALL  INIT 

C  READ  IN  ARRAY  DIMENSIONS  AND  TIME  DELAYS  FOR  THREE  SENSORS. 

101  FORMAT  ( F8 . 3 ) 

200  WRITE  (4,100) 

100  FORMAT  ( 38H  DISTANCE  BETWEEN  Ml  AND  M2  IN  METERS:) 

READ  (4,101)  R12 
WRITE  (4,102) 

102  FORMAT  (38H  DISTANCE  BETWEEN  Ml  AND  M3  IN  METERS:) 

READ  (4,101)  R13 

WRITE  (4,104) 

104  FORMAT  (34H  ANGLE  FROM  R12  TO  R13  IN  DEGREES:) 

READ  (4,101)  TH13 
WRITE  (4,106) 

106  FORMAT  (35H  ENTER  D12 ,  D13  AND  D23  IN  SECONDS:) 

READ  (4,300)  D12,  D13,  D23 

300  FORMAT  (3F8.3) 

107  FORMAT  ( F8 . 5 ) 

C  CONVERT  TH13  FROM  DEGREES  TO  RADIANS. 

TH13=TH13». 01745 

CONVERT  TIME- DELAYS  FROM  TIME  -  TO  DISTANCE  BY - - - - 

MULTIPLYING  BY  THE  SPEED  OF  SOUND. 

D12=D12*344. 

D13=D13*344. 

D23=D23*344. 

X13«R13*COS(TH13) 

Y13=R13*SIN(TH13) 

THREE  SENSOR  TARGET  LOCATION  COMPUTATION  STARTS  HERE  USING 
TIME  DELAYS  D12  AND  D13. 

A»((R12«*2)-(D12»»g))/(2.*R12) 

B=D12/R12 

A1=(X13*((D12»*2)-(R12»*2) )-R12»((D13»»2)-(R13#*2 
1)))/(2.*R12*Y13) 

B1«(D13»R12-D12*X13)/(R12»Y13) 

ZAP»SQRT( ( (2.»(A»B+A1»B1) ) »«2 ) - ( 4 . » ( B»»2+B1«»2-1 . 
1)#(A**2+A1*»2) )) 

RS1» ( -2 . * ( A*B+A1*B1 )+ZAP )/(2.*(B**2+Bl#*2-l. )) 
RS2«(-2.*(A»B+A1»B1)-ZAP)/(2.»(B»*2+B1»»2-1. ) ) 

IF(RSl.LE.O. )  GO  TO  103 
X1»A1+B1*RS1 
Y1»A+B*RS1 

THS1-(ATAN2(X1,Y1) )»57.29 


C  WRITE  OUT  DISTANCE  FROM  SENSOR  ONE  (Ml)  TO  THE 
C  TARGET  LOCATION  IN  METERS. 

WRITE  (8,105)  RSI 

105  FORMAT  ( 21H+DISTANCE  TO  TARGET  * ,  ElU.5) 

C  WRITE  OUT  ANGLE  THS1  IN  DEGREES  FROM  R12  TO  RSI. 

WRITE  (8,109)  THS1 

109  FORMAT  ( 23H+ANGLE  FROM  R12  TO  RS=,ElU.p) 

103  IF ( RS2 . LE . 0 . )  GO  TO  108 

X2=A1+B1»RS2 

Y2=A+B*RS2 

THS2  =  ( AT AN  2(X2,Y2) )*57.29 
WRITE  (8,105)  RS2 
WRITE  (8,109)  THS2 

C  THREE  SENSOR  TARGET  LOCATION  COMPUTATION  STARTS 
C  HERE  USING  TIME  DELAYS  D21  AND  D23. 

108  IF(RS1.LE.0.AND.RS2.LE.0. )  GO  TO  113 
D21=D12 

R23=SQRT((R13**2)+(R12**2)-(2.*R13#R12*C0S(TH13)) ) 
TH21=ATAN2 ( Y13 , ( R12-X13 )  ) 

X21=R12»C0S(TH2l) 

Y21*R12*SIN(TH21) 

AP=(R23**2-D23**2)/(2.*R23) 

BP=D23/R23 

A1P=(X21*( (D23**2)-(R23**2) )-R23*( ( D21*»2 ) - ( R21»»2 ) 

1)  )/(2.»R23»Y2l) 

B1P=(D21»R23-D23*X21)/(R23#Y21) 

ZAP=SQRT( ( (2.»(AP*BP+A1P»B1P)»*2)-(4.*(B»*2+ 

-  ^_^1B1P**2-1.  )*_(AP»»2+A1P**2) ) ) 

RSlP=(-2.*(AP*BP+AiP*BlP)+ZAP)/(2.*(BP**2+BlP*»2-l. ) ) 
RS2P*(-2.»(AP*BP+A1P»B1P)-Z£P)/(2.»(BP*»2+B1P**2-1.  )  ) 
IF ( RS1P . LE . 0 . )  GO  TO  110 

THS1P= ( AT AN 2 ( ( A1P+B1P»RS1P ) ,  ( AP+BP»RS1P ) ) ) *  57 . 29 
WRITE  (8,105)  RS1P 
WRITE  (8,201)  THS1P 

201  FORMAT  ( 23H+ANGLE  FROM  R23  TO  RS*,ElU.5) 

110  IF( RS2P . LE . 0 )  GO  TO  112 

THS2P=(ATAN( (A1P+B1P*RS2P) ,(AP+BP*RS2P) ) )»57 .29 
WRITE  (8,105)  RS2P 
WRITE  (8,201)  THS2P 

102  IF ( RS2P . LE . 0 . .  AND . RS1P . LE . 0 . )  GO  TO  113 
GO  TO  ilU 
113  WRITE  (8,202) 

202  FORMAT  (26H+N0  SOLUTION  WAS  OBTAINED.) 

IlU  PAUSE  10 

IF(ISENSW(17) .EQ.-l)  GO  TO  200 

STOP 

END 


C  TIME  DOMAIN  CORRELATION  PROGRAM. 

C  ROME  AIR  DEVELOPMENT  CENTER. 

C  CAPT  MIKE  CRONE  (DCTI) 

C 

C  SUBROUTINE  CORAL 

COMMON  IN ( U  ,1670 ) ,NREC ,NP 
C  IN(U, 1670-WORKING  ARRAY  WHERE  ROW  1  IS  THE 
C  FIRST  CHANNEL  OF  DATA,  ROW  2  IS  THE 

C  SECOND  CHANNEL  OF  DATA  AND  ROW 

C  3  IS  WHERE  THE  NORMALIZED  CROSS- 

C  CORRELATION  IS  STORED.  THE  ARRAY  WILL 

C  HANDLE  1670  POINTS  PER  ROW. 

C  NREC-NUMBER  OF  RECORDS. 

C  NP*NUMBER  OF  POINTS. 

B*0 

NP*NP/2 
DO  100  1*1  ,NP 
C*IN ( 1 , I ) 

100  B=B+C»*2 

DO  200  J*1 ,NP 

A1=0 

D=0 

DO  10  K*1 ,NP 
Rl-In(l.N) 

L1*K+J 

R2*IN ( 2 , Ll-1 ) 

— -  ~  ..  —  A1«R1»R2+A1  .-V  •  .•  .v; -  :  ' .  ~.t  : '  - - 

10  D»R2»*2+D 

200  IN(3,J)»A1/SQRT(B«D)*20U7.  1 

;  CALL  INIT 

;  C  CALL  INIT  INITIALIZERS  THE  CRT  DISPLAY. 

I  WRITE  (8,150) 

150  FORMAT  ( 30H  NORMALIZED  CROSS-CORRELATION*) 

|  WRITE  (8,175)  ( IN  (  3  »J )  *<7*1  ,NP ) 

1  175  FORMAT  (1H.10I6) 

|  PAUSE  1 

!  RETURN 

?  END 


